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The Higgs potential of the minimal nonminimal supersymmetric standard model (MN- 
CN \ MSSM) is investigated within the context of electroweak phase transition. We investigate 

the allowed parameter space yielding correct electroweak phase transitoin employing a high 
temperature approximation. We devote to phenomenological consequences for the Higgs 
f^*) i sector of the MNMSSM for electron-positron colliders. It is observed that a future e+e~ 

\^ \ Hnear collider with ^/s = 1000 GeV will be able to test the model with regard to electroweak 

^D . baryogenesis. 
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I. INTRODUCTION 

The baryogenesis via electroweak phase transition is one of the most interesting scenarios for 
explaining the baryon asymmetry of the universe, and it can be tested in the future accelerator 
experiments [1, 2, 3]. The mechanism for dynamically generating the baryon asymmetry at the 
electroweak scale should satisfy the three Sakharov conditions, namely, the presence of baryon 
number violating processes, the violation of both C and CP, and a departure from thermal 
equilibrium [4]. If the electroweak phase transition is responsible for the baryon asymmetry, it 
should be first order in order to satisfy the third Sakharov condition. Furthermore, the strength 
of the first order phase transition should be very strong in order to be sufficiently far away 
from thermal equilibrium. Otherwise, the baryon asymmetry might be washed out. Thus, the 
strongly first order phase transition requires that the vacuum expectation value (VEV) in the 
broken phase should be larger than the critical temperature. 

In order to satisfy the requirement that the electroweak phase transition should be strongly 
first order, the standard model is found to predict a light Higgs boson, which is much smaller 
than the experimental lower bound set by the LEP2 data [5]. On the other hand, in order 
the standard model to predict the Higgs boson mass larger than the experimental constraint 
of LEP2, the electroweak phase transition should be weakly first order. However, in this case, 
the standard model might not sufficiently generate the baryon asymmetry. Therefore, within 
the context of the baryogenesis via the electroweak phase transition, it is natural to examine 
possible extensions of the standard model to avoid this difficulty. 

One of the most plausible extensions is to embed the supersymmetry in the standard model. 
The simplest extension is called the minimal super symmetric standard model (MSSM), which 
has just two Higgs doublets in its Higgs sector. In the MSSM, the strength of the electroweak 
phase transition might be stronger than that of the standard model, at least for a specific region 
of the parameter space [6]. Thus, a strongly first order electroweak phase transition is allowed 
in this parameter region, where tan/3 is small, the right-handed stop quark is light, and one 
of the scalar Higgs boson is light. Whereas the mass of the scalar Higgs boson is above the 
experimental constraint of LEP2, the mass of the right-handed stop quark is predicted to be 
smaller than that of top quark in this parameter region. 

Introducing a Higgs singlet into the MSSM, we have the next-to-minimal supersymmetric 
standard model (NMSSM) [7], where a term which is cubic in the Higgs singlet is present in the 
tree-level Higgs potential in order to account for the breaking of the Peccei-Quinn symmetry 
which leads to a massless pseudoscalar Higgs boson. In the NMSSM, due to the presence 
of the trilinear soft supersymmetry breaking terms in the tree-level Higgs potential, the first 
order electroweak phase transition can be sufficiently strong in a wide region of the parameter 
space, where the lightest scalar Higgs boson mass can be as large as 170 GeV [8]. It is much 
easier to get a strongly first order electroweak phase transition in the NMSSM, consistent with 
phenomenological constraints, compared to the MSSM [8, 9, 10]. Also, a U(l)' extension of the 
MSSM has been investigated in Ref. [11], where a strong enough first order electroweak phase 
transition for electroweak baryogenesis is shown to exist. 

Recently a extension of the MSSM is constructed in a minimal way [12]. It is called the 
minimal nonminimal supersymmetric standard model (MNMSSM). In the MNMSSM, the cubic 
term of the Higgs singlet is absent from the superpotential of the MNMSSM. The explicit break- 
ing of the Peccei-Quinn symmetry is accomplished by the presence of the singlet tadpole term 
arising from higher order interactions. The dimensional ^ parameter of the MSSM is generated 
dynamically by means of the VEV of the neutral Higgs singlet. A number of investigations 



on the phenomenology of the MNMSSM have been done, for example, concentrating on Higgs 
production [13], spontaneous CP violation [14], and electroweak baryogenesis [15] in its Higgs 
sector. 

In this article, we consider the possibility of examining the electroweak phase transition in 
the MNMSSM at the future e"'"e~ linear colliders. First, we check if the electroweak phase 
transition in the MNMSSM may be strongly first order. We calculate the critical temperature 
where two minima of the potential are degenerate, and express the Higgs boson masses in terms 
of the critical temperature. For the parameter region where the electroweak phase transition 
is strongly first order, the production cross section for the neutral Higgs boson is calculated in 
e^e~ collisions, such as LEP2 and the future linear colliders. 

The article is organized as follows. In the next section, the Higgs sector of the MNMSSM is 
presented. In Sec. Ill, we calculate critical temperatures and the VEVs of the Higgs fields at 
them. In Sec. IV, the Higgs boson masses and their productions in e^e~ collisions are obtained 
at critical temperature. Discussions and conclusions are given in Sec. V. 

II. THE HIGGS POTENTIAL 

The Higgs sector of the MNMSSM consists of two Higgs doublets Hi (i = 1,2) and a Higgs 
singlet N, where the superpotential is given by 

W = htQH2t1i + hQHibli + XNHfeH2 , 

where the first two terms are analogous to the Yukawa interaction terms for the third generation 
quarks in the MSSM. The 2x2 antisymmetric matrix e is defined as ei2 = —£21 = 1- The third 
term corresponds to the so called /x-term of the MSSM, when the Higgs singlet develops the 
VEV. Thus, all the coupling coefficient in the above superpotential are dimensionless. There is 
no term that is cubic in the Higgs singlet in the MNMSSM. 

The tree-level Higgs potential of the MNMSSM may be decomposed into D terms, F terms, 
and soft super symmetry breaking terms as 

Vo = Vd + Vf + Vs, (1) 

where 

Vn = f{HlaH, + HlaH2f + ^f{\H2\'-\H^\^)\ 

Vs = m]j^\Hi\^ + m]i^\H2\'^ + m%\N\^-{\AxHleH2N + Yi.c.) . (2) 

Here, a are the Pauli matrices, and m^ , in'jj , and m^ are the soft SUSY breaking masses, 
and Ax is the trilinear soft SUSY breaking parameter with mass dimension. 

The above tree-level Higgs potential has a Peccei-Quinn symmetry, which leads us to an un- 
wanted axion after spontaneous symmetry breakdown, since the determinant of the pseudoscalar 
Higgs boson mass matrix is zero. In order to get rid of this axion in the MNMSSM, we introduce 
a tadpole term to the above potential due to higher-order interactions for the Higgs singlet [11]. 
The tadpole term may be given linear in the Higgs singlet as 

Utadpole = -fN - fN* , (3) 

where ^ is the tadpole coefficient. This tadpole coefficient is intrinsically a free parameter, 
but, for a phenomenological point of view, we assume that it is of the order of the soft SUSY 



breaking mass (< 1 TeV). Thus, by means of the Higgs singlet tadpole term in the MNMSSM, 
the Peccei-Quinn symmetry is explicitly broken down. 

The one-loop corrections to the tree-level Higgs potential are obtained by considering the 
effective potential method [16]. We assume that loops involving only the quarks and scalar 
quarks for the third generation are sufficient. If the scalar quark masses are degenerate, the 
one-loop effective potential due to the quark and scalar quark for the third generation is given 
as [17] 

where the soft SUSY breaking mass rh satisfies m^ = (1000 GeV)^ ^ m^ (g = t, b). In general, 
A, A\, and ^ may be complex numbers. However, as we allow no CP violation in the Higgs 
sector, they are taken real. Further, for simplicity, we take them positive, since negative values 
yield phenomenologically the same results. 

Now, in order to take into account the electroweak phase transition, we introduce the 
temperature-dependent contributions to the effective potential [18]. It is well known that in 
the standard model the high temperature approximation is consistent with the exact calculation 
of the integrals to better than 5 % for mj /T < 1.6 in the case of the fermion and for niyjjT < 2.2 
in the case of the boson [19]. Thus, we assume that the temperature for the electroweak phase 
transition in the MNMSSM is high enough as compared with the masses of the relevant particles. 
The temperature-dependent potential comes from the thermal effects due to the Higgs bosons, 
the gauge bosons, the quarks for the third generation, the lighter chargino, and three the light 
neutralinos. According to the analysis of Ref. [9], the temperature-dependent potential in the 
high temperature approximation is given by 

+ 6A2|iV|2}, (5) 

where T is the temperature of the potential and c^ = gi+ 3(72 + '^-^^- Therefore, the full effective 
potential at the one-loop level for the electroweak phase transition in the MNMSSM is 

V = Vo + ^tadpole + Vi+Vt . (6) 

which is now temperature-dependent. 

III. CRITICAL TEMPERATURE 

Now, let us evaluate the critical temperature and the VEVs of the Higgs fields at the critical 
temperature in order to examine the order and strength of the phase transition at the electroweak 
scale. We denote the temperature-dependent VEVs of the Higgs fields as vi = ((fi(T),0)) = 
{{Hf{T),H^{T))) = {H^{T)), v^ = {{0,V2{T))) = {{H+ (T), H^ {T))) = {H^iT)), and x = 
x{T) = {N(T)). Making a SU(2)i x U(l)y gauge transformation one can choose < H2 >= 
and ^2 > 0. The condition for a local minimum with < H^ >= is equivalent to the requirement 
that the charged Higgs boson have positive mass. The temperature-dependent vacuum is defined 
as the minimum of the effective potential V as 
2 2 
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Here, the soft SUSY breaking masses at the zero-temperature are given by 
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ml = -AM0)2 + ^t^(0)2A,sin2/3+^, (8) 

where tan/? = vi{Q) / V2{<d) , v'^{{)) = t;f(0) + t;i(0), and t;(0) = t;(r = 0) = 175 GeV. The 
temperature-dependent vacuum is obtained by solving three minimum equations 
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= 2mljX-2\AxviV2 + 2\^{vl + vl)x-2i^ + T'^\^x , (9) 

for the three VEVs, namely, ui, f2, x, and then substituting them into the expression for iy). 
The above minimum equations are nonlinear and thus can not be solved exactly by analytical 
method. Eventually, in (y), we are left with six free parameters as tan/3, A\, A, x(0), .^, and 
the temperature T. 

In the NMSSM, the critical temperature is defined as the largest value among the tempera- 
tures which give the second derivative of the Higgs potential zero at the origin of space of fields 
[8, 9]. In our work, we would like to use a more accurate method to determine the critical 
temperature than the method that determines it from the saddle point of the potential. Our 
method is to determine the critical temperature as the temperature at which two vacua are 
degenerate. 

Let us employ two independent methods in solving the above nonlinear equations for the 
three VEVs. First, we solve them numerically. From the third equation, one can express x at 
finite temperature in terms of the other two VEVs as 



\AxViV2 + i^ 
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Substituting this expression into the remaining two equations, and rewrite them as 

= fi{vi,V2,T) 

= f2ivi,V2,T) (11) 

where the remaining other free parameters are omitted. Solving these two equations, one can 
obtain vi and V2 for given T, and then get x. In practice, we plot /i and /2 in the {vi, f2)-plane, 
for a given temperature, by varying both vi and V2, with certain values for the remaining free 
parameters. We determine the set of points in the {vi, f2)-plane where both of /i and /2 are 
zero. In this way, we determine vi and V2 for given T, and then get x. 

In Fig. la, we plot in the (vi, t;2)-plane the contours of /i = and /2 = 0, for tan/3 = 2, 
^A = 504.5 GeV, A = 0.3, x(0) = 25 GeV, ^ = 50 GeV, and T = 100 GeV. The two contours 
intersect at five points in the (vi, V2)-pl8ine. The values of vi and V2 for these points (and the 
value of X therefrom) are the solutions of the above minimum equations. In order to examine the 
nature of these points, we plot {V) in Fig. lb for the same set of parameter values. One can see 
in Fig. lb a number of equipotential contours in the {vi, f2)-plane, and the five stationary points 
in the plane where the first derivative of (V) is zero. These points are exactly the five points 
where both /i = and /2 = 0. Among them, three points correspond to minima of the potential, 
where its second derivative of the potential is positive, and the other two points correspond to 
saddle points, where the sign of the second derivative of the potential is indeterminate. 

We compare the values of the three minima of the potential to determine the true global 
minimum, where the vacuum is defined. Among the three minima, the potential at {vi,V2,x) = 
(2.9,83.1,2.0) GeV and (426.3,385.0,234.3) GeV have the same value and lower than the one 
at the remaining point. Thus, we define two degenerate global minima, that is, the vacua, 
separated by a potential barrier. The electroweak phase transition may take place from one 
vacuum to the other through quantum tunnelling, and it is first order. The temperature which 
we set as 100 GeV is therefore the critical temperature Tc. 

Let us denote the two points in the {vi, V2, x)-space where the two degenerate vacua are de- 
fined as point A and point B. Thus, the coordinates of point A is {via, V2A,xa) = (2.9, 83.1, 2.0) 
GeV and those of point B {vib,V2b,xb) = (426.3,385.0,234.3) GeV in the {vi, V2, x)-space. We 
measure the distance between the two points in the {vi, V2, x)-space at the critical temperature 
as 

d{T,) = ^{via{Tc) - vib{Tc)Y + {v2a{Tc) - V2b{Tc)Y + {xa{Tc) - xb{Tc)Y . (12) 

As is well known, the requirement for a strongly first order electroweak phase transition is 
d{Tc) > Tc, In the present case, we have Tc = 100 GeV and d{Tc) = 569.4 GeV, which yields the 
ratio d(Tc)/Tc = 5.6. The electroweak phase transition in the present case is therefore safely a 
strongly first order one. 

We repeat the above procedure for other sets of parameter values. For simplicity, we do not 
explore the whole parameter space, but vary tan f3 and Ax while fixing the values of the other 
parameters, namely, A = 0.3, x(0) = 25 GeV, ^ = 50 GeV, and Tc = 100 GeV. Figs 2a, 3a, 4a, 
and 5a show the contours of /i = (solid curve) and /2 = (dashed curve), and Figs 2b, 3b, 
4b, and 5b show the contours of {V), for different values of tan/3 and Ax. The numerical results 
for these figures, including that of Fig. 1, are summarized by Table I. As one can see in Table 
I, the electroweak phase transitions for all the sets of parameters we consider are strongly first 
order ones. 



Now, as a second and independent way to investigate the electroweak phase transition, let 
us employ Newton's method for solving the nonlinear equations, Eq. (9). The Newton's method 
starts with a Jacobian matrix with respect to vi, V2, and x, which is defined as 
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(13) 



dxdvi ' dxdv2 ' dxdx 
The elements of the Jacobian matrix are given explicitly as 



J 



11 



2„,2 



2mjj^ + 2A^x^ + 2X^v^ + 



+ 



3hfv 



87r2 



2 + 6 log 



- 2 

m 



91+92 
2 



{3vl 



hlvj 



+ - 
rri' 



9hl 



+ ^{c^ + 3hl) 



vl 



2ht 



vi 



+ hh 



+ hlvlf 



[m' 



J' 



22 



•^33 
Jl2 
Jl3 
«/23 



2ml^ + 2A2x2 + 2\\l + ^^^(3t;i 



4„,2 



+ 



87r2 



2 + 6 log 



- 2 



2 



h't 



?;2 
^2 



+ - 






T-i2 



fc2 + 3/1?) 



2hii 



m' 



+ hj 



7)2 
^2 



[m' 



+ hhlf 



2my + 2\\vi + vi^ 

- 2X^Axx + 4A^f 1^2 

- 2A^AaW2 + 4A^t;ix 

- 2A^ylAWi + 4A^W2X 



2 I „2^ 



(ffl +52)^'1^'2 



(14) 



The Newton's method provides us with the stationary points, where the first derivative of the 
potential vanishes. The stationary points are then classified as minima, maxima, or saddle 
points, according to whether the second derivative of the potential is positive, negative, or 
indeterminate. The minima are chosen by imposing on stationary points the two conditions 
Jii > (i = 1,2,3) and det(Jjj) > 0. Among the minima, we select two points, say point 
A and point B, where {V{yiA-,V2A-,XA)) = ^{vib-,V2BiXb)) is the true global minimum of the 
potential, at a given temperature, in order to define the degenerate vacua. We check and confirm 
that the results obtained by the Newton's method coincide with the numerical results of Table 
I, for the sets of parameter values we set. 

IV. HIGGS BOSON MASS AND PRODUCTION 



In this section, we investigate the possibility of examining the parameter space where the elec- 
troweak phase transition in the MNMSSM may take place for the Higgs sector at electroweak 
baryogenesis. We calculate the neutral Higgs boson masses and their production cross sections 
in e"'"e~ collisions at the critical temperature. 

In the pseudoscalar Higgs sector, there is a neutral Goldstone boson related to Z boson. At 
the critical temperature, the elements of the 2x2 pseudoscalar mass matrix Mp{Tc) are given 
by 






2\Axx c^ 
sin 2/3 "^"6 
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Mp,,{T,;) = XvAx . (15) 

From the pseudoscalar mass matrix, the masses of the pseudoscalar Higgs bosons are obtained 
as 

1 



m\^P^{Tc 



2 



Tr Mp(Tc) T ^TiMp{TcY - 4detMp(r, 



(16) 



Note that, at zero temperature, the determinant of the pseudoscalar Higgs boson mass matrix 
is obtained as det{Mp(0)} = 2\A\^^ / sm.2(5 . If ^ / 0, the presence of the singlet tadpole term 
breaks explicitly the Peccei-Quinn symmetry. On the other hand, if ^ = 0, the singlet tadpole 
term disappears and there is a massless pseudoscalar Higgs boson at the tree level. The strength 
of the Peccei-Quinn symmetry breaking depend on the size of ^. 

For the scalar Higgs sector, the elements of the 3x3 scalar mass matrix Ms{Tc) are given 
at the critical temperature by 



Ms,,{Tc) = m'zCos'P + XAxxta.n^ + -T^ + ^T^ + ^^log 






ml 



3 j ihlml hlml 



Svr^ I fh? + m? (fh? + ttt,?)^ 



c2 2 , /i?^2 , 3/i?m? Hh^ + ml 



Ms,,{Tc) = mlsm^P + \Axxcotl3 + -T^+-^T^ + ^^^^^\og 



mt 



+ 



Ahfmt hi 



m 



Svr^ I m^ + mf [fh? + m'l)'^ 



Ms,, (T,) = ^v^Ax sin 2/3 + ^ + X^T^ , 

2 

MsATc) = {X^v^-'^)sm2p-XAxx , 
Ms^, {%) = 2X^v cos f5x - XvAx sin (3 , 

\2„ 



M523 (^c) = 2X^v sin I3x - XvAx cos (3 . (17) 

From the scalar mass matrix, the masses of the scalar Higgs bosons are obtained as 

ms.iTc? = ^Tr{M5(Te)} + 2Vl^cos|^±^| (j = 1, 2, 3) , (18) 



where 



e = cos-i' ^ 



W = -J-[Tr{M5(r,)}]2 + iTr{M5(Te)M5(r,)}, 
18 o 

U = -±-[TT{Ms{T,)}f + ^Tv{Ms{%)}Tv{Ms{T,)Ms{T,)} + 1 det{M5(r,)} . 

In e^e^ collisions, the four main production channels for the neutral Higgs bosons in the 
MNMSSM can be considered as follows [20]: 

Higgsstrahlung : e^e^ -^ ZSi , 
pair production : e~^e~ -^ SiPj , 
WW fusion : e^e~ -^ v^VeSi , 
ZZ fusion : e^e~ -^ e^e~ Si , 



where Si (i = 1, 2, 3) are three scalar Higgs bosons and Pj {j = 1, 2) are two pseudoscalar Higgs 
bosons. At the center of mass energy of LEP2, the dominant production channels for the neutral 
Higgs boson among the four channels are the Higgsstrahlung and pair production process. As 
the center of mass energy in e^e" collision increases, WW and ZZ fusion processes becomes 
important for the production of the neutral Higgs boson. We note that the Higgs boson mass 
as well as the relevant coupling for the Higgs production depends on the critical temperature. 
That is, in e^e~ collisions the production cross section for the neutral Higgs boson is evaluated 
at the critical temperature. 

We first calculate the masses of the five neutral Higgs bosons at the critical temperature. 
The results are listed in Table H, for the five sets parameters in Table I. We then calculate the 
cross sections for the production of neutral Higgs bosons via two dominant production channels 
in e~^e~ collisions with ^/s = 209 GeV (the center of mass energy of LEP2), and list the results 
in Table HI, for the five sets of parameters in Table I. We find that the 4th set of parameters 
and the 5th set of parameters yield production cross sections well below 100 fb. If we assume 
that the discovery limit of 100 fb at LEP2, the LEP2 data could not put any constraint on the 
MNMSSM with either the 4th set of parameters or with the 5th set of parameters, whereas the 
1st, the 2nd, and the 3rd set of parameters yield phenomenologically inconsistent results with 
the LEP2 data. 

We repeat the above calculation for ^ = 500 GeV (LC500) and 1000 GeV (LCIOOO). Here, 
we calculate 15 different cross sections for production of the neutral Higgs bosons via all four 
channels listed above. The results are listed respectively in Table IV and Table V. The numbers 
for the first three sets of parameters are just for comparison. We are interested in the numbers 
for the 4th and the 5th set of parameters. Assuming that the discovery limits as 20 fb for LC500, 
one can see in Table IV that, at LC500, the analysis of neutral Higgs production might determine 
whether the 4th set of parameters are acceptable experimentally. However, the neutral Higgs 
bosons might not sufficiently be produced for the 5th set of parameters, and thus the test of the 
5th set of parameters cannot be done at LC500. The MNMSSM with the 5th set of parameters 
may eventually be examined at LCIOOO, where some of the neutral Higgs boson production cross 
sections are above the discovery limit of LCIOOO, which we assume to be 5 fb. 

Now, we explore not only a number of points but a wide region in the parameter space of the 
MNMSSM whether a strongly first order electroweak phase transition is possible. The ranges 
for the parameters of the MNMSSM are set as 2 < tan/3 < 40, < A < 0.7, < ^a < 1500 
GeV, and < x(0),^ < 300 GeV. The critical temperature is determined by a temperature 
at which two degenerate global vacua exist. We search randomly in the region and discovered 
5000 points in the parameter space where the electroweak phase transition is first order by 
using the Newton's method. Among them, we obtain 1955 points where the electroweak phase 
transition satisfies d{Tc)/Tc > 1, that is, the transition is strongly first order. For each of these 
1955 points, we calculate the masses of the neutral Higgs bosons and the 15 cross sections for 
their productions at LCIOOO. The results are shown in Figs. 6 and 7. In Fig. 6, we plot the 
distributions of those 1955 points in the (r7T,5^,mp^)-plane. There is a tendency that the lower 
bound on mp-^ depends on m^j. In Fig. 7, the largest one among 15 production cross sections 
for the neutral Higgs bosons is plotted against ms^ for each of the 1955 points. For reasonable 
values of ms^, we see that at least one of the 15 production cross sections is larger than 20 fb. 
Thus, we expect that there is an ample opportunity to discover at LCIOOO at least one of the five 
neutral Higgs bosons in the MNMSSM within the context of the strongly first order electroweak 
phase transition. 



V. CONCLUSIONS 

We have investigated the possibihty of testing at the future e^e~ colhders the MNMSSM which 
accommodates the strongly first order electroweak phase transition. Both the MNMSSM and 
the NMSSM have two Higgs doublets and a singlet. But they differ in that the Peccei-Quinn 
symmetry is broken in the MNMSSM by the singlet tadpole term whereas it is broken in the 
NMSSM by a cubic term in the Higgs singlet. In the MNMSSM, we have found that 5 sets 
of parameter values for which the electroweak phase transition may take place in strongly first 
order. Radiative corrections due to the one-loop contribution from the third generation are taken 
into account, where scalar top quark masses are assumed to be degenerate. Finite temperature 
effects are considered by the one-loop finite-temperature correction in the high temperature 
approximation. Some sets of parameter values for the MNMSSM yield large enough cross sections 
for production of some neutral Higgs bosons in e+e~ collisions at y/s = 500 and 1000 GeV. By 
exploring not only the 5 sets of parameter values but a wide region in the parameter space of 
the MNMSSM, we have discovered that 1955 points in the region provide strongly first order 
transitions. For these points, we have evaluated the masses of the neutral Higgs bosons, and 
calculated the cross sections for their productions in the future e'^e~ collisions at ^/s = 1000 
GeV via four dominant channels. We have found that reasonable masses are obtained for 5*1, and 
the cross sections are considerably larger than the discovery limit of 5 fb at LCIOOO. Therefore, 
it is possible to test the MNMSSM at LCIOOO, within the context of the strongly first order 
electroweak phase transition. 
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Figure Captions 

Fig. la : Contours of /i = (solid curve) and /2 = (dashed curve) in the (fi, f2)-plane for 
tan/3 = 2 and Ax = 504.5 GeV. Other parameters are fixed as A = 0.3, x(0) = 25 GeV, and 
^ = 50 GeV. The critical temperature is set as 100 GeV. Note that the two contours intersect 
at five points in the (vi,t;2)-plane. 

Fig. lb : Equipotential contours of {V) with the same parameter values as Fig. la. There 
are five extrema; two of them are saddle points and the other three are minima. Two of the 
three minima have the same value of {V), thus define global minima or degenerate vacua. The 
coordinates in unit of GeV of the degenerate vacua at point A and point B in the (wi, t;2)-plane 
the value of x therefrom are respectively (2.9, 83.1, 2.0) and (426.3, 385.0, 234.3). 

Fig. 2a : The same as Fig. la, but with different values of tan/3 and Ax: tan/3 = 3 and 
Ax = 433 GeV. 

Fig. 2b : The same as Fig. lb, but with different values of tan/3 and Ax'- tan/3 = 3 and 
Ax = 433 GeV. The coordinates in unit of GeV of the degenerate vacua in the (wi, W2)-plane 
and the value of x therefrom are respectively (4.2, 104.8, 3.5) and (375.7, 364.9, 238.7). 

Fig. 3a : The same as Fig. la, but with different values of tan/3 and Ax', tan/3 = 10 and 
^A = 401 GeV. 

Fig. 3b : The same as Fig. lb, but with different values of tan/3 and Ax', tan/3 = 10 and 
Ax = 401 GeV. The coordinates in unit of GeV of the degenerate vacua in the (wi, W2)-plane 
and the value of x therefrom are respectively (5.0, 125.8, 10.7) and (307.7, 364.7, 360.7). 

Fig. 4a : The same as Fig. la, but with different values of tan/3 and Ax'- tan/3 = 20 and 
Ax = 499 GeV. 

Fig. 4b : The same as Fig. lb, but with different values of tan/3 and Ax'- tan/3 = 20 and 
Ax = 499 GeV. The coordinates in unit of GeV of the degenerate vacua in the (wi, W2)-plane 
and the value of a: therefrom are respectively (3.6, 128.2, 14.6) and (298.0, 409.0, 526.6). 

Fig. 5a : The same as Fig. la, but with different values of tan/3 and Ax'- tan/3 = 40 and 
Ax = 716.5 GeV. 

Fig. 5b : The same as Fig. lb, but with different values of tan/3 and Ax'- tan/3 = 40 and 
Ax = 716.5 GeV. The coordinates in unit of GeV of the degenerate vacua in the (t;i,t;2)-plane 
and the value of x therefrom are respectively (2.1, 129.0, 17.2) and (281.5, 472.7, 786.7). 
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Fig. 6 : Plot of the lightest scalar and pseudoscalar neutral Higgs boson masses for the 1955 
points in the parameter space of MNMSSM where the electroweak phase transition is strongly 
first order. The points are selected in 2 < tan/? < 40, < A < 0.7, < Ax < 1500 GeV, and 
0<x(0),C< 300 GeV. 

Fig. 7 : Plot against the mass of the lightest scalar neutral Higgs boson of the largest of the 15 
cross sections for production of any neutral Higgs bosons via any channel for the 1955 points in 
Fig. 6, in e~^e~ collisions with ^/s = 1000 GeV. 
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Table Caption 

Table I: The parameter sets of the MNMSSM where strongly first order phase transition occurs. 
Other parameters are fixed as A = 0.3, x(0) = 25 GeV, and .^ = 50 GeV. The critical temperature 
is set as 100 GeV. 



parameters 


set 1 


set 2 


set 3 


set 4 


set 5 


tan/3 


2 


3 


10 


20 


40 


Ax (GeV) 


504.5 


433.0 


401.0 


499.0 


716.5 


cor res. figures 


Figs. 1 


Figs. 2 


Figs. 3 


Figs. 4 


Figs. 5 


VIA (GeV), viB (GeV) 


2.9, 426.3 


4.2, 375.7 


5.0, 307.7 


3.6, 298.0 


2.1, 281.5 


V2A (GeV), V2B (GeV) 


83.1, 385.0 


104.8, 364.9 


125.8, 364.7 


128.2, 409.0 


129.0, 472.7 


XA (GeV), XB (GeV) 


2.0, 234.3 


3.5, 238.7 


10.7, 360.7 


14.6, 526.6 


17.2, 786.7 


d{T,) (GeV) 


569.4 


510.8 


520.6 


654.0 


887.8 


d{T,)/T, 


5.6 


5.1 


5.2 


6.5 


8.8 



Table II: The masses of the neutral Higgs bosons in the MNMSSM with the five sets of parameter 
values listed in Table I. 



mass (GeV) 


set 1 


set 2 


set 3 


set 4 


set 5 


ms^ 


46.8 


49.7 


67.6 


74.3 


76.6 


ms2 


170.6 


165.3 


158.2 


156.7 


156.3 


ms,, 


297.8 


251.6 


220.8 


297.0 


478.1 


"^Pl 


65.0 


61.6 


70.3 


75.2 


77.0 


mp2 


298.8 


252.9 


221.3 


297.0 


478.0 
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Table III: Cross sections for production of the neutral Higgs bosons in e+e^ collision with ^/s = 
209 GeV via the Higgsstrahlung and the pair production channel with the five sets of parameter 
values listed in Table I. 



cross section (fb) 


set 1 


set 2 


set 3 


set 4 


set 5 


<y{zSi) 


718.2 


674.4 


264.2 


79.8 


20.0 


(yiSiPi) 


161.9 


193.5 


34.0 


2.9 


0.1 



Table IV: Cross sections for production of all neutral Higgs bosons in e+e~ collision with ^/s 
= 500 GeV via the Higgsstrahlung, the pair production, the WW fusion, and the ZZ fusion 
channel, with the five sets of parameter values listed in Table I. 



cross section (fb) 


set 1 


set 2 


set 3 


set 4 


set 5 


<y{ZSi) 


53.5 


51.0 


22.7 


7.3 


1.8 


o{ZS2) 


5.6 


4.0 


0.03 


0.01 


0.01 


o^ZS^) 


1.7 


4.6 


24.8 


19.4 


0.0 


o^S^Px) 


30.8 


36.4 


8.7 


0.9 


0.06 


o^S^Px) 


25.0 


15.1 


0.3 


0.005 


0.0001 


<y{S^Pi) 


0.5 


2.7 


9.0 


1.8 


0.0 


o^SxPi) 


0.9 


3.1 


8.4 


1.7 


0.0 


a{S2P2) 


0.1 


0.7 


0.2 


0.004 


0.0 


ct{SsP2) 


0.0 


0.0 


3.6 


0.0 


0.0 


a{UeJ^eSl) 


115.6 


108.4 


43.3 


13.4 


3.4 


a{9eVeS2) 


5.4 


4.0 


0.03 


0.02 


0.01 


<j{i?eVeSz) 


0.6 


2.5 


17.1 


7.7 


0.002 


a{e+e-Si) 


12.3 


11.6 


4.6 


1.4 


0.3 


a{e+e~S2) 


0.5 


0.4 


0.003 


0.002 


0.001 


<j{e+e-S^) 


0.07 


0.2 


1.8 


0.8 


0.0002 
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Table V: The same as Table IV, but with a higher center of mass energy, y/s = 1000 GeV. 



cross section (fb) 


set 1 


set 2 


set 3 


set 4 


set 5 


(t{ZSi) 


10.7 


10.3 


4.6 


1.5 


0.3 


a{ZS2) 


1.4 


1.0 


0.008 


0.004 


0.003 


cr{ZSs) 


0.7 


1.6 


7.5 


9.1 


6.2 


criSiPi) 


7.7 


9.1 


2.2 


0.2 


0.01 


a{S2Pi) 


8.3 


4.9 


0.09 


0.001 


4.5 


criS^Pi) 


0.4 


1.4 


3.8 


1.4 


0.2 


<y(SiP2) 


0.7 


1.5 


3.6 


1.4 


0.2 


ct{S2P2) 


0.7 


0.8 


0.1 


0.01 


0.0006 


(t{S3P2) 


0.03 


0.2 


6.0 


8.1 


0.4 


a{Uel^eSl) 


238.3 


225.5 


95.7 


30.3 


7.7 


a{i?eJ^eS2) 


19.4 


14.2 


0.1 


0.06 


0.04 


a{VeT^eSz) 


6.9 


16.5 


85.9 


79.5 


30.0 


a{e+e-Si) 


26.3 


24.9 


10.6 


3.3 


0.8 


a{e+e~S2) 


2.1 


1.6 


0.01 


0.007 


0.005 


a{e+e-S^) 


0.7 


1.8 


9.7 


9.1 


3.4 
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Fig. la: Contours of /i = (solid curve) and /2 = (dashed curve) in the (fi,f2)-plane for 
tan/3 = 2 and A\ = 504.5 GeV. Other parameters are fixed as A = 0.3, x(0) = 25 GeV, and 
^ = 50 GeV. The critical temperature is set as 100 GeV. Note that the two contours intersect 
at five points in the (fi,t;2)-plane. 
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Fig. lb: Equipotential contours of {V) with the same parameter values as Fig. la. There 
are five extrema; two of them are saddle points and the other three are minima. Two of the 
three minima have the same value of {V), thus define global minima or degenerate vacua. The 
coordinates in unit of GeV of the degenerate vacua at point A and point B in the (vi, t;2)-plane 
the value of x therefrom are respectively (2.9, 83.1, 2.0) and (426.3, 385.0, 234.3). 
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Fig. 2a: The same as Fig. la, but with different values of tan/3 and A\: tan/3 = 3 and 
Ax = 433 GeV. 
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Fig. 2b: The same as Fig. lb, but with different values of tan/3 and Ax: tan/3 = 3 and 
A\ = 433 GeV. The coordinates in unit of GeV of the degenerate vacua in the (t;i,t;2)-plane 
and the value of x therefrom are respectively (4.2, 104.8, 3.5) and (375.7, 364.9, 238.7). 
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Fig. 3a: The same as Fig. la, but with different values of tan/3 and A\: tan/3 = 10 and 
Ax = 401 GeV. 
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Fig. 3b: The same as Fig. lb, but with different values of tan/3 and Ax: tan/3 = 10 and 
A\ = 401 GeV. The coordinates in unit of GeV of the degenerate vacua in the (wi, W2)-plane 
and the value of x therefrom are respectively (5.0, 125.8, 10.7) and (307.7, 364.7, 360.7). 
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Fig. 4a: The same as Fig. la, but with different values of tan/3 and Ax: tan/3 = 20 and 
Ax = 499 GeV. 
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Fig. 4b: The same as Fig. lb, but with different values of tan/3 and Ax: tan/3 = 20 and 
A\ = 499 GeV. The coordinates in unit of GeV of the degenerate vacua in the (t;i,t;2)-plane 
and the value of x therefrom are respectively (3.6, 128.2, 14.6) and (298.0, 409.0, 526.6). 
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Fig. 5a: The same as Fig. la, but with different values of tan/3 and A\: tan/3 = 40 and 
Ax = 716.5 GeV. 
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Fig. 5b: The same as Fig. lb, but with different values of tan/3 and Ax: tan/3 = 40 and 
A\ = 716.5 GeV. The coordinates in unit of GeV of the degenerate vacua in the (t;i,t;2)-plane 
and the value of x therefrom are respectively (2.1, 129.0, 17.2) and (281.5, 472.7, 786.7). 
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Fig. 6: Plot of the lightest scalar and pseudoscalar neutral Higgs boson masses for the 1955 
points in the parameter space of MNMSSM where the electroweak phase transition is strongly 
first order. The points are selected in 2 < tan/3 < 40, < A < 0.7, < Ax < 1500 GeV, and 
0<x(0),C< 300 GeV. 
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Fig. 7: Plot against the mass of the Ughtest scalar neutral Higgs boson of the largest of the 15 
cross sections for production of any neutral Higgs bosons via any channel for the 1955 points in 
Fig. 6, in e'^e" collisions with y^ = 1000 GeV. 
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